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1 . Introduction 

One of the important object of DUMAND is the measurement of high 
energy neutrinos. And also in the underground laboratories, such as 
Gran Sasso, neutrino physics will become more and more important stfcject. 
In such a case atmospheric neutrinos behave sometimes as guest and 
sometimes become an obstacle. In this paper focussing on high 

energy neutrinos 1 TeV) , we carried out a new calculation of 
atmospheric neutrino intensities taking into account EMI effects observed 
in P-A collisions by accelerator, recent measurement of primary cosmic 
ray spectrum and results of cosmic ray muon spectrum and charge ratio. 
Another features of the present calculation are a) taking into account 
kinematics of three body decays of kaons and charm particles in diffusion 
equations and b) taking into account energy dependence of kaon production. 


2. Methods of the calculation 

A . 

Energy moment in which particle h collide with nucleus A and 
produce particle c is given as following; 
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On the other hand in nucleus -nucleus collisions, we take A-A collision 

as sum of "nucleon''-A collision. Where "nucleon" is considered that 

a part of nucleon is quark- gluon (Q-G) state. So Z ^ is modified 

pc 

taking into account the primary cosmic ray composition as follows; 

Z'^ =6 Z ^ + fl-6 )[fl-e) Z ^ + 6 Z ^ 1 . 

pc p pc ^ p-^ ■' pc ^ pC -• ’ 

where 6^ is the proton excess at the top of the atmosphere and the value 
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is 0.79. 3 is the rate of Q-G state and 0.237 estimated from cosmic 

ray muon charge ratio. More details of the above formulas and the 

A ll 

values of can be found in Minorikawa and Mitsui ^ . 

Next we consider prompt neutrino production. Prompt neutrinos are 

urn “ 4 * •• 

produced through decays of charm particles, such as pp'^ DDX or A^DX. 
Inclusive distribution of D-meson is very similar to the distribution of 
mesons produced by h-A collision and the distribution of baryon is 
similar to that of proton without diffractive part. From the above 
points, energy moments of charm particles are presented as following; 
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where 3 = a <n >A /a <n >A . As a cross section, decay mode and 

branching ratio of charm particles, we take those presented by 
21 

Castagnoli et al ^ and the collision mean free path of charm particle 

2 2 + 
with nucleon was taken 200 g/cm for D-meson and 100 g/cm for A^ 

baryon at the energy of 1 TeV. For neutrino intensities produced by 

three body decays of kaon and charm particles we can calculate using 

diffusion equation taking into accotmt kinematics by Hagedom ’ and 

given as follows; 
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where is decay constant and = 850, Bj^o=201.8, B^o=9. 04x10 , 

7 8 ^ 

Bjp± = 4.33x10 and B^+ = 2.12x10 in GeV. Focussing on the neutrino 
from the decay D ->■ Kp§, we present as following; 
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Primary spectrum was taken as 1.87 E~ (E^^ 5 x 10 GeV/nucleon) and 
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191 (Eq> 5 X lO^GeV/nucleon) in unit o£ (an^sec sr GeV)~^. 

5. Results 

Calculated differential neutrino spectra are presented in Fig.l as 
solid lines. The calculations performed in the early, Volkova^^ and 
Inazawa and Kobayakawa^^ , are shown in the figure , however latter 
intensity was raised by a factor it for kinematical correction. 

By integrating differential neutrino spectrum from 0° to 360° for 
azimuthal angle and from 0° to 90° for zenith angle, downward going 
integral neutrino intensities were estimated and presented in Fig. 2. 

4. Discussion 

Rate of the intensity of electron neutrino to that of muon neutrino 
is 5.5% at 1 TeV and 2.6% over 10 GeV. And also rate of the intensity 
of neutrino to that of anti -neutrino is 2.2 for muon neutrino and 1.7 for 
electron neutrino, independently of energy. The integral intensity of 
the prompt neutrino exceeds the conventional one at about 10 TeV for 
electron neutrino, and in the energy region of several hundreds TeV for 
muon neutrino. Because the above energy region will be covered by 
DUMAND, it seems to be able to detect the prompt effects. 

This calculation was performed using the computer FACOM M180-II AD of 
the Institute for Nuclear Study, University of Tokyo. 
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Fig. 1 : 

Differential neutrino 
spectrum. 

0° and 90° mean zenith angle. 



Fig. 2 : 

Neutrino intensity I (>E) 
obtained by integrating the 
differential spectrum from 
0° to 360° for azimuthal angle 
and from 0° to 90° for zenith 
angle. 









